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Abstract 

We study the transverse quantum ANNNI model in the region of high frustration (k > 0.5) 
using the DMRG algorithm. We obtain a precise determination of the phase diagram, showing 
clear evidence for the existence of a floating phase, separated from the paramagnetic modulated 
phase by a high-order critical line ending at the multicritical point. We obtain simple and accurate 
formulae for the two critical lines. 

PACS numbers: 75.10.Jm, 73.43.Nq, 05.10.Cc 
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I. THE MODEL 



The ANNNI model is an axial Ising model with competing next-nearest-neighbor antifer- 
romagnetic coupling in one direction. It is a paradigm for the study of competition between 
magnetic ordering, frustration and thermal disordering effects. 

In the Hamiltonian limit, we consider a one-dimensional quantum spin 5* = | chain 
interacting with an external magnetic field, called the TAM model (transverse ANNNI). 

The TAM Hamiltonian for L spins with open boundary conditions reads^'^ 

L-l L-2 L 

i=l i=l i=l 

We use the "traditional" notation k — —J^l^x- The notations A = JxjB and V — B are 
sometimes used in the literature. 

The sign of J\ is immaterial, since the Hamiltonian is invariant under the transformation 

<^(-i)vf. (2) 

Likewise, the sign of B is immaterial. Without loss of generality, we set J\ — \. We restrict 
ourselves to positive k and even L. 

We also consider fixed boundary conditions, where we add to the extremities of the chain 
two fixed spins (Tq and (Tl+i, with the possibilities of parallel (cro,cri,+i =tt) or antiparallel 
(aojCL+i =ti) boundary conditions. 

In the region of high frustration > 0.5), despite extensive studies'^"^, the phase diagram 
of the transverse ANNNI model is not well known. For low the model is known to be in 
the gapless "antiphase" ttii- It undergoes a second-order phase transition at a magnetic 
field Bi{k,). The existence of a "fioating" phase, massless and with slowly decaying spin 
correlation functions, up to a Kosterliz-Thouless phase transition at a magnetic field B2{k), 
is an open question. For high B, the TAM is known to be in a paramagnetic modulated 
phase. 

II. OBSERVABLES 

We measure the two lowest energies Eq and Ei, the mass gap A — Ei — Eq, the entangle- 
ment entropy Sa (see below), and two spin-spin correlation functions: the "slow" correlation 
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function 

Csid) = <Ji/2+i+,), l<d<L/2 (3) 

and the "fast" correlation function^^ 

c/(^^) = K/2-d<^!/2+i+d), 0<d<L/2. (4) 

Interesting quantities related to the correlation functions are: 
the overlap o of Cs{d) with the antiphase correlation function 

2 

cM = (-1)L(<^-V2)/2J^ 0= - J]c,(ci)c„(d); (5) 

d=l 



the average fast correlation function (times an oscillating sign) 

cf = {-ir^'j—^Y.^fidy, (6) 

the range of the fast correlation function 



A. Entanglement entropy 



It is possible to study an order-disorder phase transition using the entanglement entropy^. 
We divide the system of size L into a left subsystem of size i and a right subsystem of size 
L — £, and define 

SA{e;L)^-Tr{pAlnpA), (8) 

where A denotes the degrees of freedom of the left subsystem, B the degrees of freedom of 
the right subsystem, and pa — T^b |^o)(^o|; note that SA{i; L) — Sa{L—£; L). For a critical 
system we expect (neglecting lattice artifacts) 

Sa{£] L) - (c/6) log(Lsin(7r£/L)), (9) 

where c is the conformal anomaly number (central charge) of the corresponding conformal 
field theory, and ~ means "up to a (non- universal) additive constant" ; for the case of interest 
for the infinite- volume DMRG, £ = |L, and sin(7r£/L) only shifts the additive constant. For 
a noncritical system, we expect 

Sa{\L-L) - (c/6)[logL + s(L/e)], (10) 



where s{x) is a universal finite-size scaling function satisfying the constraints s(0) = and 
s{x) ~ — \ogx for large x. 



B. Domain- wall energy 

So far, we only considered open boundary conditions. Following Rcf. 9, we define the 
domain-wall energy (note that our definition of L differs by 2 from the definition of Ref. 9) 



where Eq and Eq are the ground state energies with parallel and antiparallel boundary 
conditions respectively. 

III. THE ALGORITHM 

We implement the density matrix renormalization group (DMRG) algorithm described in 
Ref. 10. We sample the Ug lowest energy levels with equal weights, i.e., we use the reduced 
density matrix 



(see Eq. (26) of Ref. 10). Usually, since we are interested in the mass gap A, we set Us = 2. 
We identify system and environment (for antiparallel boundary conditions, up to a spin fiip 
erf — >■ —df, (T? ~crf)- The typical dimensions of the truncated system and environment 
M range from 80 to 160; in the following, M — 80 will be understood, unless M is explicitly 
quoted. 

The crucial part of the numerical computation is finding the lowest eigenvalues and 
eigenvectors of the superblock Hamiltonian; we employ the Implicitly Restarted Arnoldi Al- 
gorithm implemented in Arpack^^, in the routine dsaupd used in mode 1. We use (typically) 
100 Lanczos vectors and require convergence to machine precision, obtaining residual norms 



We observe a truncated weight (the sum of the eigenvalues of the density matrix whose 
eigenvectors are dropped in the truncation) e ~ 10~^ for "normal" configurations, and 
e ~ 10~^ for peaks of A (see below). 



^Dw(«:,i^,^) = (-l)''/'+' eI\k,B,L)-EI^{k,B,L) , 



(11) 




(12) 



Hx-Xx\/\X 
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We managed to diagonalize the system exactly up to L = 22; both the finite- and the 
infinite- volume DMRG algorithm reproduce the results of exact diagonahzation. For moder- 
ate L, finite- and the infinite- volume DMRG give consistent results. For higher L, discrepan- 
cies between finite- and infinite- volume DMRG and M-dependence of A becomes noticeable; 
they are strongly observable-dependent, and they will be discussed below, where results on 
observables are presented. 

During a run of the finite- volume DMRG algorithm on a system with L„ sites, information 
about the system/environment for all smaller system is available. It is therefore possible, 
with a moderate extra numerical effort, to estimate the observables for all the systems 
with L < Ln sites. These estimates almost coincide with the results obtained running 
independently at each L. The additional errors introduced by this procedure will also be 
discussed below. 

The finite-volume DMRG algorithm at large L requires a very large amount of memory; 
however, since the observables at each lattice size are accessed only twice per cycle, they can 
be conveniently kept on disk, requiring only a very large amount of disk space; for L — 600 
and M — 80, e.g., ca. 6 Gbytes are required. 

IV. PHASES AT K = 0.75 

We will first focus our attention on the model at k — 0.75, and later on extend the study 
to other values of k. 

Running the infinite-volume DMRG algorithm at k, — 0.75 and B < 0.257, with open 
boundary conditions, we observe that the mass gap A vanishes exponentially in L, apart 
from numerical errors due to the fact that A is computed as Ei — Eq: sec, e.g., Fig. 1. 

The slow correlation function Cs{d) almost coincides with Ca{d); the overlap o approaches 
a value very close to 1 with corrections proportional to 1/L. The fast correlation function 
Cf{d) is constant and close to ±1, apart from d — and ds close to L/2; the range R is 
almost exactly L/4 (the value for a constant Cf{d)) and the average Cf approaches a value 
very close to —1 with corrections proportional to 1/L. 

There is a very sharp phase transition at 0.257 < Bi < 0.258. We will postpone its 
detailed study, since it is best done using £^dw- 

Running at k = 0.75 and B > 0.258, the mass gap A as a function of L at fixed B shows 



5 




sharp peaks with a frequency increasing with B: see Fig. 2. Each peak match exactly a 
change of sign of c/. At each S, for L smaller than the first peak we observe signals very 
similar to the case B < 0.257; for higher L, we observe that the minima of A seem to go 
to zero for B < 0.4 and to a nonzero limit for B > 0.5; however, the determination of A 
from the infinite-volume DMRG is not accurate for L > 200; finite-volume DMRG data with 
M — 80 become unreliable for L > 300; we show in Fig. 3 the case B — 0.3. 

We performed a finite-size analysis of A, similar to the analysis of Ref. 2, but with a 
complication arising from the peak structure. We run the finite-volume DMRG algorithm 
for B = 0.4, 0.41, 0.42, 0.43, 0.44, 0.45, 0.46 and L„ > 292 corresponding to a minimum 
of A. For each B, we select the minima of A and define Al{B, k) outside the minima by 
interpolation in L. We now take two values Li and L2 and look for the intersection Bi{Li, L2) 
of the two curves LiAl^ {B, k) and I/2 A^^^ {B, k) vs. B (interpolating in B at fixed L and k as 
needed). The results are shown in Fig. 4: we note that M — 120 and M — 160 data almost 
coincide, and even M — 80 data are adequate in the range of Ls considered; we quote as a 
final result B2 — 0.424(3). The data presented here were obtained from a run at a single L„ 
for each B (see Sect. Ill); in order to check that the error introduced is under control, we 
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FIG. 2: Mass gap peaks in the L-B plane for k = 0.75. 
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FIG. 3: Minima of mass gap A vs. L for n = 0.75 and B = 0.3. 



also performed separate runs for all the values of L required, for = 80 and for = 120 
at i? = 0.42, and repeated the analysis: the values of Bi{Li, L2) never change by more than 
0.0005. 
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FIG. 4: The intersection Bi vs. L, for k = 0.75. 

For 0.258 < B < 0.45, the slow correlation function Cs{d) at fixed L shows oscillations 
with power-law damping, in rough agreement with 

Cs{d) = ad-'' cos{qd + 0), (13) 

with r/ ~ for i? = 0.258, increasing with B but remaining smaller than |. For B > 0.5, 
Cs{d) at fixed L shows oscillations with exponential damping. We tried to extract rj by fitting 
c^(d), smoothed by taking a running average over [27r/g + |J points, to the form ad~'^^. In 
Fig. 5 we show the typical case B = 0.425; c1{d) for different values of L converge not to 
a single curve but to two separate curves, with different values of t], preventing a precise 
determination of rj. 

The range of the fast correlation function R should distinguish clearly the floating phase, 
where i? ^ oo as L ^ oo (since 77 < |) , from the paramagnetic phase, where R has a flnite 
limit as L 00. A flrst problem is the presence of oscillations, with dips corresponding to 
the peaks of A (see Fig. 6); it is solved by selecting the values of R at the Ls corresponding 
to the peaks of A. After this operation, R vs. L at flxed B and M is well fltted to the form 

U(L) = ; (14) 

if we plot the asymptotic slope sr vs. B we should be able to see a drop towards in cor- 
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data, L = 252, 266, 280, 294, 308, 322, 336, 350, 364, 378, 392, 406 

data, L = 418, 432, 446, 460, 474, 488, 502, 516, 530, 544, 558, 572 
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FIG. 5: The smoothed squared slow correlation function (?g{d) for k = 0.75 and B = 0.425, from 
the finite-volume DMRG at L„ = 710 and M = 120, for values of L corresponding to minima of 
A. (Data at M = 80 give a very similar plot.) 




9 
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FIG. 7: The entanglement entropy Sa{\L; L) vs. the size of the system L, for k = 0.75, B = 0.425, 
Ln = 710, and infinite-volume DMRG. 

respondence with the phase transition. However, in the critical region, R (and Sr) strongly 
depend on M: see Fig. 6; we can only conclude that B2 — 0.45. 

The analysis of the range of the slow correlation function gives similar, but even less 
precise, results. We conclude that the spin correlation functions are unsuitable for the 
precise determination of the floating/paramagnetic phase transition. 

A. Entanglement entropy 

We show in Fig. 7 the entanglement entropy Sa for the typical case B = 0.425. Finite- 
volume DMRG essentially reproduces the results of infinite-volume DMRG at the same M. 
We can estimate that the DMRG determination of 5*^ is reliable up to L = 200 for M = 80, 
up to L = 300 for M = 120, and up to L = 400 for M = 160. For each value of M, within 
the given range of L, the difference between finite- and infinite-volume DMRG results is less 
than 0.001. We can therefore compute Sa using the faster infinite-volume DMRG, and this 
allows us to work at larger values of M. 

In the antiphase, Sa is essentially constant, indicating a very small correlation length. 
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FIG. 8: The reciprocal correlation length 1/^ (determined from Sa) vs. B, for k = 0.75. Fitting 
the M = 120 and M = 160 data for 20 < L < 200 we obtain results almost identical to the 
M = 80 results plotted; likewise, fitting the M = 160 data for 20 < L < 300 we obtain results 
almost identical to the M = 120 results plotted; we omit these results from the plot for readability. 
All data for B < 0.43 are consistent with zero, with a very small error 
at the scale of the plot. 



~ 10 ®) which is not visible 



The simple Ansatz for the finite-size scaling function entering Eq. (10), 

s(x) = -ln(a; + e-"^) (15) 

with a = 1, is found to fit the entanglement entropy data very well (excluding just the very 
smallest lattices with L < 10) in all cases for the floating and paramagnetic phases; c is 
always compatible with 1. The best determination of ^, obtained by fitting Sa with c = 1 
fixed, is shown in Fig. 8; our final estimate is B2 = 0.44(1). 



B. Domain- wall energy 

So far, we only considered open boundary conditions. We now switch to fixed boundary 
conditions, in order to compute the do main- wall energy -Edw! with fixed boundary condi- 
tions, there are no problems with quasi-degenerate energy levels (typically, A > 0.01) or 
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FIG. 9: The infinite-volume domain-wall energy E^o vs. i?, for k = 0.75. 

peaks in A associated with level crossings, and truncated weights are e ~ 10^^ or smaller 
for M = 80. Even if we are not interested in the mass gap, we run with Ug = 2, which gives 
results more stable than rig = 1. 
We may fit -Edw to the form 



Ebw = aexp{—dL)L " + 



in the antiphase and 



E, 



DW 



a exp{—dL) 



cos{kL + 0)1 — I sin(A;L + 0)| 



(16) 



(17) 



in the floating phase (with d = 0) and in the paramagnetic phase^. The fits, excluding 
(typically) lattices with L < 16, are of very good quality and stable. 

Eq. (16), fits perfectly E^w for B < 0.257, giving E^o — > 0, constant u ^ 1.6 and 
d ~ 0.008 for B / Bi. Eq. (17), with d = 0, fits perfectly Edw for B > 0.257 and gives 
/c — i> for B \ Bi. The best estimators of Bi are E^o in the antiphase and k"^ in the floating 
phase, both vanishing linearly at Bi, see Figs. 9 and 10; the final estimate of the critical 
field is Bi = 0.2574(2). 

So far, we obtained results very similar to those of Ref. 9. We turn now to the problem 
of identifying the floating phase, i.e., a region with d = 0. The data generated with the 
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FIG. 10: The squared modulation parameter fc^ vs. B, for k = 0.75. 



infinite-volume DMRG at 0.3 < B < 0.4 seem to indicate c/ < 0; this appears to be an 
artifact of the infinite-volume DMRG, as we can see from the comparison of ii^Dw evaluated 
with the finite- and infinite-volume DMRG at B = 0.3, shown in Fig. 11. 

We must therefore resort to the resource-consuming finite- volume DMRG. We checked in 
several instances that M = 80 is sufficient to obtain accurate results and that obtaining the 
data for all Ls from a run at a single L„ is acceptable: in all cases, the determinations of d 
are well within the error quoted. 

We show the results in Fig. 12. The value of d obtained from the finite- volume DMRG is 
consistent with zero up to i? = 0.425, where we estimate ^ = l/d> 10^. The smoothness of 
d vs. B suggest a higher-order, possibly Kosterliz-Thouless, phase transition. Note that d is 
quite compatible with 1/^ of Fig. 8 (apart from a normalization stemming from the different 
definition of correlation length), and so is the resulting B2 = 0.435(10). 

The precise determination of the transition point of a Kosterliz-Thouless phase transition 
is a notoriously difficult problem: it is always possible that a system with a huge correlation 
length is mistaken for a critical system^^. Indeed, the correlation length is expected to 
diverge very rapidly when the Kosterliz-Thouless critical coupling is approached; following 
Ref. 13, we could conjecture a behavior ^ oc exp(&/ {B — B2)) for B \ B2. A fit to the above 
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FIG. 11: Determinations of the domain-wall energy -EdW) multiplied by L, vs. L, for k = 0.75 and 
B = 0.3: finite-volume DMRG at M = 80, 120, 160; infinite-volume DMRG at L„ = 368, M = 80; 
fit to Eq. (17) of the Ln = 368 data. 
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FIG. 12: The decay parameter d of the domain-wall energy -EdW) computed with the finite-volume 
DMRG, vs. B, for k = 0.75. 
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form of the data of Fig. 12 gives unstable results for B2, indicating that the determination 
of B2 from -Bdw should be treated with some caution. On the other hand, a fit to the data 
of Fig. 8 gives results which are stable and fully consistent with the estimate of B2 from Sa 
given above. 

The two determinations of B2 from and Sa, obtained by quite different methods, are 
accurate and in agreement with each other (and with the less reliable determination from 
E'dw); moreover, for each method, we see no trend of B2 decreasing with increasing L or 
M (see especially Figs. 4 and 8). We can conclude that, while it is possible that the errors 
on B2 are underestimated, it is difficult to believe that the floating phase might disappear 
completely on larger systems. 

V. PHASE DIAGRAM 

The study of the phase transitions at other values of ks is very similar to the one at 
K, — 0.75 presented in Sect. IV and we can avoid repeating the details. We selected for our 
analysis the values k = 0.5, 0.52, 0.55, 0.6, 0.75, 1.0, 1.25, 1.5, 2.0, 5.0. 

At K = 0.5, the DMRG algorithm becomes inefficient at low B, and we are unable to 
run at B < 0.01. We see no sign of a floating phase: the curves LAl{B,k,) vs. B almost 
coincide for B < 0.06, and the intersections are very unstable. Determining ^ by fitting Sa 
with c = 1 fixed, we see no sign of ^ = 00, see Fig. 13; given the poor convergence in M for 
small B, we estimate B2 < 0.04. The analysis of -Edw does not give precise results for B2. 
The behavior of the modulation parameter k (see Fig. 14), which goes to a nonzero value as 
B ^ 0, hints at the very peculiar nature of the multicritical point at, k = 0.5, B = 

For K, = 0.52 and 0.55, the quality of the determinations of Bi and of the determination 
of B2 from Sa is similar to those at n = 0.75; on the other hand, the analysis of and 
£^DW do not give precise results for B2. We present the plot of 1/^ tor k — 0.52 in Fig. 15: 
the difference from Fig. 13 is remarkable. 

For 0.6 < K, < 1.5, there are no relevant differences from the case k, — 0.75 described in 
Sect. IV; we only present the determinations of Bi and B2 in Table I. For k — 2, the only 
difference is that -Edw is not fitted well by Eq. (17) in the floating and paramagnetic phases, 
and therefore the determination of B2 from -Edw is unreliable. 

In the case k — 5, the determination of Bi and B2 is rather imprecise: E'dw is not fitted 
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FIG. 13: The reciprocal correlation length 1/^ (determined from Sa) vs. B, for n = 0.5. 
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FIG. 14: The squared modulation parameter k'^ vs. B, for k = 0.5. 



well by Eq. (16) in the antiphase, and it is fitted poorly by Eq. (17) in the floating and 
paramagnetic phases; it is very hard to get precise results from Al, since the modulation 
parameter is very small {k < 0.01 in the floating phase). It is still possible to estimate Bi 
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FIG. 15: The reciprocal correlation length 1/^ (determined from Sa) vs. B, for k = 0.52. 



K 


Bi (-Bdw) 


Bi (o) 


B2 (Al) 


B2 {Sa) 


B2 (Ebw) 


0.5 








< 0.06 


< 0.04 


< 0.08 


0.52 


0.0201(1) 


0.021(1) 


0.095(15) 


0.115(5) 


0.12(3) 


0.55 


0.0501(2) 


0.052(2) 


0.160(15) 


0.175(5) 


0.18(2) 


0.6 


0.1015(2) 


0.103(2) 


0.235(6) 


0.25(1) 


0.25(1) 


0.75 


0.2574(2) 


0.2575(5) 


0.424(3) 


0.44(1) 


0.425(10) 


1.0 


0.5213(2) 


0.522(2) 


0.700(5) 


0.72(1) 


0.71(1) 


1.25 


0.7867(2) 


0.785(2) 


0.972(4) 


1.00(1) 


0.98(1) 


1.5 


1.0514(2) 


1.045(5) 


1.235(3) 


1.26(1) 


1.26(1) 


2.0 


1.5775(2) 


1.576(2) 


1.756(6) 


1.79(1) 


1.79(3) 


5.0 




4.667(3) 




4.88(1) 





TABLE I: Determinations of the transition fields Bi and B2 by different techniques. 

from o and B2 from Sa- 

For all the values of k, considered, the different determinations of Bi and B2 are in 
substantial agreement with each other: this is a strong argument supporting the reliability 
of our results. It should be noticed, however, that the determination of B2 from A^, is 
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FIG. 16: Phase diagram in the k-B plane; the solid lines for k > 0.5 correspond to Eq. (18). 

systematically lower than the determination from Sa, possibly indicating that the error on 
the determination from reported in Table I is underestimated. 

We can beautifully summarize all the above results by noticing that all the determinations 
of Bi and B2 are consistent with 

5i(k) = 1.05(k- i), B2{k) = 1.05^{k - ^){k - 0.1). (18) 

Finally, we draw the phase diagram in the k-B plane in Fig. 16. The region k < 0.5 was 
studied in Ref. 2; the critical line separating the paramagnetic modulated and paramagnetic 
unmodulated phases is known analytically^^. The data in the region k > 0.5 are taken from 
the present work; note that earlier results^"^ provided only a qualitative picture of the phase 
diagram in this region. 

A very interesting question is whether the floating phase extends up to k = 00 or it 
terminates at finite k; we found that the floating phase extends at least up to k = 5. 

VI. SUMMARY AND CONCLUSIONS 

We applied the DMRG algorithm to the study of the quantum transverse ANNNI model 
in the region of high frustration [k > 0.5). 
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We obtained clear evidence for the existence of a floating phase for n > 0.5, extending at 
least up to K = 5. The floating phase is separated from the paramagnetic modulated phase 
by a high-order (possibly Kosterliz-Thouless) critical line, ending at the multicritical point 
{k = 0.5, -8 = 0); the corresponding central charge is c = 1. In Ref. 7, the floating phase 
was shown to have a finite extent at k = 0.5; our study cannot exclude a floating phase of 
very small extent, i.e., < B2{n^0.5) < 0.04. 

We obtained precise estimates for the critical points, verifying that different methods give 
consistent results. Simple and accurate formulae for the two critical lines are reported in 
Eq. (18). 

Very helpful discussions with Walter Selke and Pasquale Calabrese are gratefully acknowl- 
edged. 
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